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^ ■ The propagation of surface acoustic waves in a rotating anisotropic 

crystal is studied. The crystal is monoclinic and cut along a plane 
'^ I containing the normal to the symmetry plane; this normal is also the 

^ ' axis of rotation. The secular equation is obtained explicitly using the 

Q ■ "method of the polarization vector", and it shows that the wave is 

dispersive and decelerates with increasing rotation rate. The case of 
orthorhombic symmetry is also treated. The surface wave speed is 

^ • computed for 12 monoclinic and 8 rhombic crystals, and for a large 

range of the rotation rate/wave frequency ratio. 
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1 Introduction 

Introduced more than thirty years ago, Surface Acoustic Wave (SAW) de- 
vices have been used with great success by the telecommunication industry: 
nowadays, they are produced in large quantities (several billions per year) 
and used in wireless transmission and reception technology for color televi- 
sion sets, cell phones, global positioning systems, etc. In recent years, new 
applications for SAW devices have emerged, namely acoustic sensors, which 
are passive (no power supply is needed), resistant, almost non-aging, cheap 
(only one photo-lithographic process is involved in the production), light (less 
than Ig) and can be operated remotely and wirelessly. For instance [Ij SAW 
identification tags are used for highway toll collection in Norway and for the 
Munich subway system; SAW temperature sensors can achieve a resolution 
of 0.02°C from — 196°C up to 500°C; wide ranges and fine resolutions are also 
achieved for pressure, torque, or current sensors; etc. Also, the automotive 
industry is engaged in the search for an "intelligent tire" which could provide 
direct information on its current state as the car is moving; in this context 
SAW sensors have been used to measure tire pressure [2] or friction [3j as the 
wheel rotates. In general, SAW devices may be used as angular rate sensors 
(gyroscopes) to measure frequency shifts due to the rotation [U El |6] . In the 
present paper, an investigation of the effect of rotation upon the speed of 
surface (Rayleigh) waves in an anisotropic crystal is presented. 

The crystal may possess as little as a single plane of symmetry. It is cut 
along any plane containing the normal to the symmetry plane and is assumed 
to rotate at a constant rate about this normal. The surface wave is polarized 
in the symmetry plane. In other words, it suffices to consider the propagation 
of a surface wave in the Xi direction of a monoclinic crystal with symmetry 
plane at X3 = 0, cut along the X2 = plane, and rotating about the 0:3- 
axis (see Fig. 1). The secular equation for rotating materials was obtained 
by others but in simpler settings: by Clarke and Burdess in an isotropic 
material, first for small rotation rate/wave frequency ratios [1], then for any 
ratio [5]; by Grigor'evskii, Gulyaev, and Kozlov [7] also for isotropic materials 
but neglecting the centrifugal force; and by Fang, Yang, and Jiang [6J for 
crystals having tetragonal symmetry. Here, the analysis is fully developed 
for crystals with a single symmetry plane, up to the derivation of the secular 
equation in explicit form, that is an equation giving the Rayleigh wave speed 
in terms of the elastic parameters and of the rotation rate. 

The equation is reached in Section 3, after the governing equations have 
been written down in Section 2. The secular equation turns out to be a 
polynomial of degree 8 for the squared wave speed and also for the squared 
rotation rate/wave frequency ratio. In the simpler case of orthorhombic sym- 



metry (Section 4), the polynomial is of degree 6. The Rayleigh wave speed 
is computed numerically for 20 specific anisotropic materials (12 monoclinic, 
8 orthorhombic) and for a rotation rate/wave frequency ratio varying from 
to 10. Of course, this range is way beyond the elastic behavior limit, and is 
irrealistic for pratical purposes where the frequency of a SAW device is typi- 
cally in the lOOkHz-lOMHz range. It is presented to show that the method of 
resolution is exact and not approximate, applies for any rate of rotation, and 
that in contrast with the non-rotating case, the secular equation is dispersive. 
At small rotation rates, and for certain crystals such as PZT-5, other papers 
[HI [8] show that the Rayleigh wave speed may at first increase slightly with 
the rotation frequency /wave frequency ratio. At large ratios, it is seen here 
that the wave speed decreases with increasing ratios. These variations are 
crucial to the understanding and correct design of rotating SAW sensors or 
SAW signal processing devices. A recent article [9] describes the manufac- 
turing of a 1 cm X 1 cm SAW gyroscope and how the rotation rate may be 
measured using SAW technology. Another example that springs to mind is 
that of "spinning missiles" [TU] for which it is reasonnable to speculate that 
the communication is ensured via SAW generation and processing of high- 
frequency signals modified by the rotation. Finally in Section 5, the merits of 
several methods of derivation for the secular equation in non-rotating crys- 
tals are discussed. This paper aims to provide a theoretical and analytical 
framework for the study of surface acoustic waves in rotating crystals. 

2 Basic Equations 

We consider a half-space 0:2 > occupied by a homogeneous anisotropic crys- 
tal possessing one plane of symmetry at 0:3 = 0, and rotating at a constant 
angular velocity fl about the xa-axis. We study the propagation of a surface 
(Rayleigh) wave in the Xi-direction, with attenuation in the a;2-direction. In 
the rotating Cartesian frame {Oxi, 0x2,0x3) = (0,i,j,k), the equations of 
motion are |TT] 

div (7 = putt + 2pfik X u.i + pi7^k x (k x u), (1) 

where cr is the Cauchy stress tensor, p is the constant mass density of the 
material, and the comma denotes differentiation. The second term in the 
right hand-side of Eq.([T]) is due to the Coriolis acceleration, the third is due to 
the centrifugal acceleration. Note that Eq.([T]) represents the time- dependent 
part of the full equations of motion. The time- independent part, namely div 
cr"* = pfl'^kx [kx (u'^+x)], where u* = u''(x) and af^- = Cijkiulj^ must be solved 
separately. The questions remain of (a) whether an actual time-independent 



solution exists for all Q and if it does, of (b) whether the boundary conditions 
of a traction-free rotating half-space may be satisfied without perturbating 
the time-dependent boundary value problem. These questions do not seem 
to have been addressed in the literature, but some preliminary work seem 
however to suggest that (a) and (b) may be answered positively, at least 
within the framework of small amplitude waves superimposed upon a large 
elastic deformation. 

Now, turning back to the time-dependent equations ([T]), the mechanical 
displacement u is taken in the form 



u{xi,X2,X3,t) = \J{kx2)e 



ik{xi—vt) 



(2) 



showing a sinusoidal propagation with speed v and wave number k in the Xi- 
direction, and the possibility of an attenuation in the a;2-direction through 
the unknown function \J{kx2)- 

We wish to describe the influence of the frame rotation upon the speed 
of Rayleigh waves, and to this end, we introduce the following quantities. 



X 



pv 



n/{kv) = n/co, 



where a; is the real frequency of the wave. 

For two-dimensional motions {du/dx^ = 0) such as Eq.([2]), the anisotropy 
of a crystal possessing X3 = as a symmetry plane is described by the 
following strain-stress relationship |12j : 
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where the strain components e^j are defined in terms of the displacement 
components by: 2ejj = Uij + Uj^i, and the s^ • are the reduced compliances. 



Alternatively, the equivalent strain-stress relations can be used [121 P-39]) 



C°e°, C's' = 1, 



(3) 



where 



[crii,(T22, 0-23, 0-31,(112] 



[eii,...,2ei2] 
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The C°- are ele- 



ments, in the Voigt notation, of the fourth-order elastic stiffness tensor Cijki- 
Table 1 shows the relevant reduced compliances of 12 different monoclinic 
crystals, computed from the corresponding stiffnesses as collected by Chad- 
wick and Wilson |13|; the last column gives the corresponding Rayleigh wave 
speed in the non- rotating case [H]. 



In view of the form Eq.(|5]) for the displacement, we introduce the functions 
ti, ^2 for the tractions cri2, cr22 on the planes X2=const. as, 

(Ti2(xi,X2,X3,t) = ikti{kx2)e'''^''^~'"''\ (T22{xi,X2,X3,t) = ikt2{kx2)e'^^''^'''^\ 

Then, substituting Eqs.(|2]) and ^ into the equations of motion ([T]), we derive 
the following system of linear first-order differential equations for [/i, f/25 ^i, 

*2, 



U' 

t' 



N, 



N3 + 

T 



(4) 



2j , t = [^1,^2]^! and the prime denotes differentiation with 



Ni 

'l + 5^)Xl Ni 

where \J = [U 

respect to kx2- Note that, as in the static case [15], the anti-plane strain 
(stress) decouples from the plane strain (stress) and need not be considered 
for this problem. This decoupling would not occur if the the crystal was 
rotating about the xi-axis or the X2-axis 0. 

The surface X2 = is free of tractions and so, the boundary conditions 
are 

ti(0)=t2(0) = 0. (5) 

In Eq.(j4]), Ni and N2 are the same as the 2x2 submatrices of the 6x6 
fundamental elasticity matrix N from Ingebrigsten and Tonning [TB]. Their 
real matrix N3 however has been modified by the introduction of off-diagonal 
pure imaginary terms. Explicitly, we have 
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where the quantities r2, 
parameters as [TH |17] 



rg, n22, ^26, n^Q, rj are given in terms of the elastic 



1 
'11 



nm 



'11 



re 



'11 

'16 



'16 



'11 



'16 



^2 



"-22 



'11 



I 
12 

/ ' 
11 

■^12 



■^12 
■^22 



^^26 



'11 



■^12 



'16 
'26 



Thus the rotation of the crystal perturbs the equations of motion in three 
ways: the introduction of dispersion through 5; a shift of magnitude S^ in 
X = pv"^ for the lower left submatrix of N proportional to the 2x2 unit 
matrix; and the modification of N3, which is diagonal in the non-rotating 
case (note that the new matrix N3 is Hermitian: N3 = N3 .) Despite these 
modifications, the secular equation can be obtained explicitly for the surface 
wave speed, using a method proposed by Currie [TB| and by Taziev [12] 
for non-rotating anisotropic crystals with and without a plane of symmetry, 
respectively. 



3 Secular Equation 

The method of the polarization vector was first presented by Currie [18] to 
derive the secular equation for Rayleigh waves in the symmetry plane of 
monoclinic non-rotating crystals. Then Taziev [19] generalized the method 
to triclinic (no symmetry plane) crystals. This method takes advantage of 
the Cayley-Hamilton theorem for the fundamental matrix N, which implies 
that only n matrices N'^ {k = 1, . . . ,n) are linearly independent (n = 3 for 
monoclinic crystals, n = 5 for triclinic crystals). Currie used the matrices 
N, N^, N^; Taziev, the matrices N, N^, N^, N^ N^ Recently, Ting [20] 
placed their results within the context of the Stroh-Barnett-Lothe formal- 
ism and improved on them by showing that the choices of N~^, N, N^ for 
monoclinic crystals and of N~^, N~^, N, N^, N^ for triclinic crystals lead 
to simpler and more explicit secular equations. His approach is now adapted 
to our present context of a rotating crystal with one symmetry plane. An 
alternative derivation, not based on the Stroh-Barnett-Lothe formalism, is 
available elsewhere [211 [22] • 

We seek solutions to the equations of motion Eq. (jl]) presenting exponen- 
tial decay with distance 

V{kx2) = ae''=^"^ tikx2) = be''=P"^ S(p) > 0, 

where the constant vectors a and b are related through [121 p. 139]: 6j = 
{Ckiii + pCi2k2)(^k- Then the equations of motion Eq.(jl]) give 



P 



N 



(6) 



where N is the 4x4 matrix in Eq.(jl]). This eigenvalue problem yields a 
quartic for p. We limit our investigation to the subsonic range, defined as 
the greatest interval of values for v where the determinant of N— pi possesses 
two roots pi, p2, with positive imaginary parts. We call ai, a^, and bi, b2, 
the vectors a and b corresponding to each root. Then the solution is of the 
form pSl p. 141] 

U = A < e^^P* > q, t = B < e'''P* > q, < e^^^* >= diag (e^'^Pi^^ e*^^^^^), 

where A = [ai,a2], B = [bi,b2], and q is a constant vector. Using the 
boundary conditions Eq.©, we have at the free surface X2 = 0, 

Bq = 0, and u(xi,0,X3,t) = aRe^'^^^i-^*^ aR = Aq. (7) 

Moreover, the matrices A and B satisfy the orthogonality condition 

— X — T 

B A + A B = 0. 
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Now, the eigenrelation Eq.OH]) may be generalized for any positive or neg- 
ative integer n as 



P" 



N" 



where N" 



n(") 



N 



(n) 



K(") Ni' 



(n)T 



(say). 



(9) 



Exphcitly, the elements of N" are computed by multiplication of N or its 
inverse by itself. For instance, K*^") for n = 1,2, —1 is given by K^-*-) = 
N3 + X(l + (52)1, 
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where Z? is a real denominator common to the K^, whose expression is too 



(-1) 



long to reproduce and which turns out to be irrelevant for the derivation of 
the secular equation. 

Now we write in turn the second vector line of Eq.([n])i for pi and for p2, 
and deduce 

K(")A + n/"^B = B diag ipi,p2). 

Multiplying this equahty to the left by slr^ = q^A and to the right by q, 
and using Eqs.([7j),([8]), we conclude that (see [20j for the non-rotating case). 



At n = —1, 1, 2, and hr = [1, a]^ (say), three equations follow: 
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We re- arrange this system as: Fi^gk = hi,hj introducing the following quan- 
tities, 

Fu = D^{k[-''^), Fu = D'=s{Ki-''^), F,s = DKi-'\ 

F21 = 0, F22 = s[,%kS), F23 = s',,kS, 

F31 = s;iK(E:g)), Fs2 = s[,^iKil^), F33 = 0, 

gi = a + a, g2 = i{a — a) , 5'3 = cta, 

h = -DK[^'\ h2 = -s'nK'i'^, h = -s'nK[f. 

Note that the explicit expressions for the non-dimensional quantities Fj^ and 
hi in terms of X = pv"^, 6 = fl/u, and the s[j are easily read off Eqs. flT(I]) . flTT]) . 
For instance, F12 = 2s[2SX, F32 = —2s'iqSX, h2 = I — s[i{l + S'^)X, and so 
on. 

The linear non-homogeneous system Fg = h has a unique solution for g. 
Introducing A = det F and A^ {k = 1,2,3), the determinant of the matrix 
obtained from F by replacing its k-th column with h, we write the solution 
as gk = Afc/A. But the components of g are related one to another through 
5'3 = (fl'i/2)^ + ((72/2)^. This relation is the explicit secular equation for 
Rayleigh waves on an anisotropic crystal rotating in its plane of symmetry, 

A^ + A^ - 4A3A = 0. (13) 

This equation is a polynomial of degree 8 in X = pv"^, and also of degree 8 
in 6^. Because 6 = fl/u appears only in even powers in the secular equation, 
the Rayleigh speed obtained as a root of Eq.f lT^ does not depend on the 
sense of rotation. Numerically, we find that the rotation slows the Rayleigh 
wave down and that the speed is a monotone decreasing function of 6. We see 
this behavior on Fig. 2, where the dependence of the Rayleigh wave speed 
upon 6 is shown for the 12 monoclinic crystals from Table 1. The curves 
are arranged in the same order as in the Table, from the slowest (diphenyl, 
starting at 1276 m/s) to the fastest (diallage, starting at 4000 m/s). 

The secular equation is valid for any crystal possessing at least one plane 
of symmetry, as long as the half-space is cut along a plane containing the 
normal to the plane of symmetry. In particular, it is also valid for orthorhom- 
bic crystals when the plane of cut contains one of the crystallographic axes. 
When this plane contains two crystallographic axes, the secular equation 
factorizes and a separate treatment is required. 



4 Orthorhombic Materials 

When the material possesses three orthogonal planes of symmetry and the 
axes (0,i,j,k) are aligned with the crystallographic axes, some compliances 
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vanish: s[q = s'2q = 0. Table 2 lists the values of the relevant reduced 
compliances for 8 rhombic crystals, computed from the corresponding values 
of the stiffnesses as collected by Shutilov [21]. The corresponding Rayleigh 
wave speed vr in the non-rotating case (last column) is found from the exact 
secular equation [2] 



(1 - Xs[,)y/1-Xs',,- X^s[,[s',, - X{s[,s',, - s\%)] = 0. 
When the frame is rotating, the system Fg = h reduces to 
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i^isl 
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\h{] 
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-^23 
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h2 
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where F31 7^ and 

F12 = 2s[26X, Fi3 = 1 - {s[, + 46)(1 + 6')X + s[,s',,{l - 6'YX^, 

F22 = -2s'n6X, F23 = s\,{l + 6^)X, 

h = [^'22(1 + 5') - (s'n42 - 4)(1 - S'fX]X, h2 = l- s[,il + 6')X. 

From this new system of equations, we deduce that (?2 = ^2/^ and g^ 
A3/A, where 



FToFo-i — FooF^ 



22-^13, 



A, 



^1-^23 — ^2-^13, ^3 — -^12^2 " -^22^1, 



and also that gi = a + a = Q, implying that (72 = 2m, (73 = —a^ = ((72/2)^ 
as well. This last equality is the explicit secular equation for Rayleigh waves 
on an orthorhombic crystal rotating in one plane of symmetry, 

Aj - 4A3A = 0. 

This equation is a polynomial of degree 6 in X = pv"^ and in 6 = Q/u. As 
in the monoclinic case above, the roots are even functions of S. Numerically, 
the results are similar to those of the monoclinic case, as Fig. 3 shows for 
the eight orthorhombic crystals of Table 2. Again, the curves are arranged 
in the same order as in the Table, from the slowest (sulfur, starting at 1628 
m/s) to the fastest (benzophenone, starting at 4723 m/s). 



5 Concluding Remarks 

Several methods have been proposed to derive explicitly the secular equa- 
tion for surface waves in non-rotating monoclinic crystals with the plane 
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of symmetry at 0:3 = 0. This Author [14J wrote the equations of motion 
as a system of two second-order differential equations for the tractions t: 
(iiktk ~ iPikt'k ~ liktk = 0, whcrc o., (3, 7 are 2x2 real symmetric matrices. 
Then the method of first integrals [25] yields the secular equation. The equa- 
tions of motion Eq. (j4j) may also be written in a similar manner for a rotating 
crystal, but a, (3, 7 become complex and the method of first integrals is 
no longer applicable as such. Next, Ting [TTj assumed an exponential form 
for t{kx2) and obtained the secular equation through some simple algebraic 
manipulations, taking advantage of the fact that ai2 = /322 = 0; in the ro- 
tating case however, these quantities are no longer zero. Furs [26] (using the 
displacement field) and this Author [27] (using the traction field) devised 
yet another method, where the secular equation is the resultant of two poly- 
nomials; again, having real quantities for the components of N is a crucial 
property, no longer true for rotating crystals. 

All in all, it seems that the method of the polarization vector is the most 
appropriate for the case of a rotating crystal. Note that a simple derivation 
of its main result (IT^ . not relying on the Stroh formalism, was presented 
recently 
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Figure 1: Monoclinic crystal with symmetry plane at X3 = 0, cut alon^ 
X2 = 0, and rotating about x^ at constant angular velocity Q. 
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Table 1. Values of the reduced compliances (10 ^^ m^/N), density (kg/m^) 
and (non-rotating) surface wave speed (m/s) for 12 monoclinic crystals. 



material 


^ii 


■^22 


S'l2 


■^16 


•^26 


*66 


P 


VR 


diphenyl 


854 


1858 


-366 


-698 


-1.44 


5049 


1114 


1276 


tin fluoride 


345 


228 


-59.2 


-197 


120 


922 


4875 


1339 


tartaric acid 


343 


211 


-164 


-223 


301 


1650 


1760 


1756 


oligoclase 


133 


227 


-108 


97.0 


-160 


483 


2638 


2413 


micro cline 


94.5 


165 


-35.1 


47.2 


1.69 


446 


2561 


2816 


gypsum 


243 


130 


-68.6 


32.9 


28.1 


326 


2310 


3011 


hornblende 


63.3 


103 


-32.7 


-15.8 


-2.72 


320 


3120 


3049 


aegirite-augite 


53.6 


78.4 


-21.0 


-10.6 


-33.5 


237 


3420 


3382 


epidote 


53.3 


49.6 


-11.3 


17.7 


-3.74 


237 


3400 


3409 


augite 


54.5 


64.4 


-19.5 


-19.0 


-15.7 


211 


3320 


3615 


diopside 


53.1 


58.6 


-20.1 


24.0 


6.98 


186 


3310 


3799 


diallage 


49.8 


69.1 


-11.3 


-6.88 


-14.5 


166 


3300 


4000 
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Figure 2: Rayleigh wave speeds for 12 monoclinic crystals rotating about 0:3. 
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Table 2. Values of the reduced compliances (10 ^^ m^/N), density (kg/m^) 
and (non-rotating) surface wave speed (m/s) for 8 orthorhomhic crystals. 



material 


^ii 


■^22 


^'12 


■^66 


P 


vr 


sulfur 


65.1 


76.2 


-42.2 


132 


2070 


1628 


iodic acid 


36.1 


20.1 


-7.88 


57.5 


4630 


1678 


a-uranium 


4.89 


5.29 


-1.13 


13.5 


19000 


1819 


rochelle salt 


49.3 


33.0 


-18.2 


102 


1775 


2114 


sodium-tartrate 


32.1 


27.1 


-16.8 


102 


1818 


2197 


strotium formate 


24.5 


30.9 


-7.32 


58.1 


2250 


2451 


olivine 


3.26 


5.34 


-0.97 


12.6 


3324 


4599 


benzophenone 


13.0 


13.9 


-7.17 


27.9 


1219 


4723 
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Figure 3: Rayleigh wave speeds for 8 rhombic crystals rotating about x^. 
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